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ABSTRACT 

We present a list of "best possible" estimates of low-degree p-mode frequencies, from 8640 
days of observations made by the Birmingham Solar-Oscillations Network (BiSON). This is 
the longest stretch of helioseismic data ever used for this purpose, giving exquisite precision 
in the estimated frequencies. Every effort has been made in the analysis to ensure that the 
frequency estimates are also accurate. In addition to presenting the raw best-fitting frequencies 
from our "peak bagging" analysis, we also provide tables of corrected frequencies pertinent 
to the quiet-Sun and an intermediate level of solar activity. 
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1 INTRODUCTION 

It has been known for over forty years that the su rface of the Sun 
oscillates with a period of around 5 minutes (see iLeighton et al.l 
1 19621 : lEvans & Micharoj[l962l) . Originally thought to be local in 
nature, observations of the 5-minute oscillations using unimaged 
sunlight (i.e., Sun-as-a-star observations) proved t hat the oscilla- 
tions were actually a global phenomenon (see e.g., IClaverie et all 
1 19791) . The oscillations are the detectable manifestation of acous- 
tic waves that are generated by turbulent motion in the solar con- 
vection zone. Waves of certain frequencies interfere constructively 
leading to a large number of resonant p modes, so called du e to the 
gradient of pressure being the dominant restoring force fsee lUlrichl 
ll97CtlLeibacher & Steinll97lh . 

In Helioseismology, estimates of the frequencies of solar 
p modes are used as inputs to the inverse problem of trying 
to better constrain the par ameters of solar models (see e.g., 
IChristensen-Dals gaard 2002]). Hence it is important to determine 
the frequency estimates to high precision and high accuracy. 

The different resonant modes are characterised by three inte- 
gers: the radial order, n, angular degree, t, and azimuthal order, 
m. Low-£ modes travel deep into the solar interior, while higher- 
degree modes are constrained within the outer regions. Hence, by 
determining frequencies of a full range of modes it is possible to use 
inversions to determine the structure of the solar interior as a func- 
tion of depth. However, the constraints on these inferences become 
poorer the deeper one goes into the interior as there are fewer and 
fewer modes available which penetrate to successive depths. Hence 
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when dealing with low-^" modes, which are observed most effec- 
tively by global (i.e., Sun-as-a-star) observations, obtaining precise 
frequency estimates is particularly desirable. 

The precision with which any p-mode frequency can be de- 
termined is directly related to the length of observations one 
makes. The Birmi ngham Solar-Oscillations Network (BiSON) (see 
IChaplinet all 19961) has collected more than 30 years of observa- 
tions of global helioseismology data and is therefore in a position 
to provide estimates of low-/' p modes to unprecedented levels of 
precision. In Section [2] we give details of the BiSON network and 
explain how the BiSON time series are generated. We also give the 
parameters of the particular data set used in obtaining the p-mode 
frequencies given in this paper. 

The most common method of determining the frequencies, 
which is used here, is to perform a Fourier transform of the time se- 
ries in order to generate the frequency power spectrum. The modes 
are then parameterized by fitting an asymmetric Lorentzian model 
to the peaks, the parameters being frequency, height, width, rota- 
tional splitting and fractional asymmetry. This fitting is colloquially 
referred to as "peak bagging" analysis. 

In order to give accurate estimates of these parameters one 
needs a model that matches the underlying limit spectrum of the 
mode peaks as closely as possibl e. Here we have employed the 
"pseudo-global" fitting model of iFletcher et all fe008t l2009al lbT). 
which we describe briefly in Sectionf3] 

The one drawback of using very long time series in order to 
determine highly precise frequency estimates is the presence of the 
solar cycle. During one activity cycle of the Sun, the p-mode fre- 
quencies change by up to 1 /jHz (many times the precision we are 
attempting to report the frequency estimates to). However, it is pos- 
sible to track the change in frequency over time and therefore the 
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solar cycle effect can be accounted and corrected for. We present 
raw fitted frequencies and frequencies corrected for the solar cycle 
in Section|4](for the cases of the quiet-Sun and intermediate levels 
of solar activity). 



2 THE BISON DATA 

The BiSON network is comprised of six semi- and fully-automated 
solar observing stations that are dedicated to collection of low- 
degree (Sun-as-a-star) helioseismic data. The stations are situated 
at various sites around the world in order to provide as continuous 
observations as possible. The six sites are at: Mount Wilson in Cal- 
ifornia, Las Campanas in Chile, Izana in Tenerife, Sutherland in 
South Africa, Carnarvon in Western Australia and Narrabri in New 
South Wales. 

At each of the six stations a resonance scattering spectrome- 
ter (RSS) (see Brook es~et al.lll978h is used to measure the Doppler 
velocity shift of the 770-nm Dl potassium absorption line. This is 
done by passing the incident solar light through a cell containing a 
vapour of potassium atoms at a temperature of about 100°C. Pho- 
tons of the appropriate energy are resonantly scattered by atoms in 
the cell. As the light from this process is emitted isotropically it is 
possible, by placing detectors at right angles to the incident beam, 
to record only those photons which have undergone resonance scat- 
tering. The absorption cross section of atoms in the vapour cell is 
much narrower than that of the solar Fraunhofer line because the 
temperature is lower and there is no rotational broadening. There- 
fore the intensity of the recorded light will be proportional to the 
intensity of light emitted from a narrow band of the solar absorption 
line. 

The intrinsic sensitivity of the RSS is improved by moving the 
passbands out onto the wings of the Fraunhofer line where the slope 
is greatest and hence any given line shift will give a greater change 
in measured intensity. The vapour cell is placed in a longitudinal 
magnetic field, which causes Zeeman splitting to occur. The single 
line is thus split into a multiplet with a separation dependent on the 
field strength. The splitting alters the atomic state in such a way that 
atoms will interact with circularly polarized light. The blue-shifted 
transition is sensitive to one hand of circular polarization whilst the 
red-shifted transition is sensitive to the other. 

By using a suitable combination of a linear polariser and a 
quarter-wave plate the incident light can be circularly polarised and 
hence, by switching quickly between one-handedness and the other, 
it is possible to measure the light intensity in the blue wing, I b and 
red wing, I r , almost simultaneously. From these measurements a 
ratio, R, is formed which gives a near-linear proxy for the velocity 
shift of the solar line: 



Finally, to obtain velocity measurements of the solar oscilla- 
tions, a third-order polynomial function of the station-Sun line-of 
sight velocity is fitted to R. The oscillations signal is then recovered 
by subtracting R from the polynomial function and calibrated using 
the fitted gradient of R versus station velocity. 

These measurements of the solar oscillations are collected at 
each station and combined to form a continuous time series with 
a duty cycle of around 86% (for 2007 and 2008). Izana has been 
collecting data since 1976, hence there are over 30 years of BiSON 
data to call upon. However, with only one station collecting data 
and at only certain times of the year the fill for these early observa- 
tions is only of the order of 10%. During the mid 1980's and early 



1990's the remaining BiSON station's were gradually brought on- 
line and the window function improved to around 80% and above. 

In this paper we have chosen to work with an 8640-day set of 
40-s cadence BiSON data running from 4th April 1985 until 16th 
December 2008. The start date corresponds to a time just after the 
third BiSON station came into use in Carnarvon, Western Australia. 
This means the data set starts with a fill of around 20-30% over the 
first few years which steadily improves to around 80% by 1994 
once all six stations were full operational. 



3 THE "PSEUDO-GLOBAL" TECHNIQUE 

The best way of determining estimates of the p-mode frequencies 
is to take the Fourier transform of the time series in order to get the 
frequency power spectrum. For Sun-as-star observations the power 
spectrum has a very distinct pattern. A series of sharp peaks is seen 
that coincide with the various p-mode frequencies, with the largest 
peaks signifying the strongest modes at around 3000//Hz (i.e., a 
period of around 5 minutes). 

On closer inspection the peaks are seen to be grouped in pairs. 
These pairs consist alternately of t = and 2 modes and i = 1 and 
3 modes. In power spectra of long BiSON time series it is also pos- 
sible to make out the strongest £ = 4 modes, which tend to have 
frequencies a little lower than the 1 = 2 modes (thus the t = 0/2 
pair could be considered an t = 0/2/4 triplet). Even closer inspec- 
tion reveals that that all but the I = modes are rotationally split. 
However, even with a very high resolution spectrum this splitting 
cannot be seen at high frequencies since modes in this regime have 
short lifetimes and consequently their peaks have large widths in 
the power spectrum. 

To determine the frequency estimates an asymmetric 
Lorentzian model is fitted to the various peaks in the frequency 
power spectrum. The traditional way of doing this has been to di- 
vide the frequency power spectrum into a series of fitting regions 
and fit the modes in pairs (i.e., I = and 2 mode s together and 
C = 1 and 3 together) (see e.g., IChaplin et al.ll 19991) . This enables 
one to determine how the mode parameters depend on both fre- 
quency (overtone number) and angular degree without the need to 
fit the entire spectrum simultaneously. 

However, recent work has shown that t his method can result 
in systematic bias in the retu rned frequencies jjimenez-Reves et"al] 
l2008l : iFletcher et alj |2009a). The reason for this is that the model 
only accounts for modes whose central frequencies are within the 
fitting region. Therefore, power from modes whose central frequen- 
cies lie outside this region will not be accounted for. This imperfect 
match between the fitting model and the underlying profile can lead 
to bias in the mode parameters. 

A global fitting approach where one fits the entire spectrum 
simultaneously would overcome this problem. However it would 
mean using a fitting model with many hundreds of parameters 
which, when fitting a frequency power spectrum made up of more 
than 18 million points, leads to prohibitively long computing times. 
Also, such a complicated model may lead to an increased possibil- 
ity of premature convergence. 

There fore, in order t o elimin ate the bias the "pseudo-global" 
method of Fletc her et al. (2009a) was employed. This approach 
works by first setting up a model for the full spectrum (FSM) using 
the parameter estimates returned from the traditional "pair-by-pair" 
method. Then, as with the pair-by-pair approach, the frequency 
power spectrum is divided into a series of fitting regions. Each re- 
gion is then fitted using the FSM, but only the parameters of the 
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modes within each particular fitting region are allowed to vary. In 
this way we limit the number of parameters being varied at any one 
time, while still including information in the model from all the 
modes. A much more detailed account of this fitting method can be 
found in lFletcher etail J2009ah . 

Two further improvements were made in the pseudo-global 
model used here compared with the techniques used when deter- 
mining previously published BiSON frequencies. The first is the 
inclusion of t = 4 modes in the model. Secondly we take account 
of the gaps in the BiSON data by convolving the fitting model with 
the spectral window, a s opposed to simply fi tting for the prominent 
diurnal sidebands (see lFletcher et al. 2009b). 

When employing the pseudo-global method the fitting regions 
were 130yuHz wide and centered at the midpoint of the £ = 0/2 and 
£ = 1/3 pairs. The C = 4 modes were also included in the model 
when these modes were visible. Within each fitting region the fol- 
lowing equation is used to fit the power spectral density, P, as a 
function of frequency, v: 



P(V): 



where, 



P'(v) + n(v) + ^ 

tin 



f (m A,{\ + 2bx em ) 



* w. 



2[v - (vQf + mSvf)] 



(2) 



(3) 



W is the spectral window function and the asterisk denotes a con- 
volution. P' is the power from the remainder of the power spectrum 
outside of the fitting region. A ; denotes the mode heights and f lm the 
relative m component fractional height ratios (assumed to take fixed 
values), b is the fractional asymmetry of the modes (fixed across the 
fitting region), vOf are the mode frequencies, Sv/ are the rotational 
splittings and T( are the mode widths (fixed a cross m compo nents). 
Finally, n is the background term given by the Harvey ( 1985) model 
to account for the solar granulation with a constant term, /? to ac- 
count for instrumental noise. 



n(v) 



2cr 2 T 



1 + (2tcvt)< 



+ P, 



(4) 



where t and a represent the time constant and standard devia- 
tion respectively of the granulation signal. Within the fitting region 
only the /? term is allowed to vary. The values of cr and t were 
fixed beforehand, by fitting the Harvey model to the low and high- 
frequency regions of the frequency power spectrum (i.e., outside 
the region of the modes). 

Once all the fitting regions had been covered, an improved 
set of parameter values (compared to those returned from the tradi- 
tional fitting method) will have been determined for the entire spec- 
trum. These new values can then be used to set up another model 
for the entire spectrum and the whole pseudo-global process can be 
run again. In general the process only needs to be repeated one or 
two times before no further significant improvements are seen in 
the mode parameters between successive iterations. 

Differences between the frequency estimates from the pseudo- 
global technique and the traditional approach are small. However, 
because of the high precision of the BiSON frequencies the differ- 
ences are significant and also systematic. The frequencies from the 
pseudo-global technique are between 0.02 - 0.04 /jHz (1-2 sigma) 
lower than those determined from a traditional approach over the 
region 2200 - 3500 fiHz. 



4 A DEFINITIVE LIST OF BISON FREQUENCIES 

The raw, best-fitting frequencies from our pseudo-global peak- 
bagging analysis are listed in Table 1. We have augmented this list 
at very low frequencies with previously published BiSON frequen- 
cies that come fro m work where we de l iber ately targeted the low- 
frequency regime ( Chapl in et al. I I2002U2007I) . We have included 
one additional frequency to those previously published, the n = 7, 
C = 1 at 1185.592 yiiHz, however only the m = -1 component was 
detected. To estimate the frequency centroid we therefore added to 
the frequency of the observed component the mean rotational fre- 
quency splitting of nearby non-radial modes. 

Tables of frequencies have to be used with care, because p- 
mode frequencies vary with the level of solar activity. We therefore 
provide sufficient information to enable our frequencies to be cor- 
rected to an activity level of choice. This will help comparisons and 
combinations with observations made by other instruments, which 
do not span the 8640-day period covered here. We present two ta- 
bles of frequencies, which were produced by applying suitable cor- 
rection procedures to t he raw frequencies listed in Table 1. We use 
the 10.7-cm radio flux (Tapping & DeTracev 1990) as our adopted 
proxy of the global level of solar activity. We have not corrected the 
i = 4 frequencies, for the reasons described near the end of this sec- 
tion. Corrections have not been applied to the very low-frequency 
modes (<1450 ^/Hz), as these show very little cycle variation and 
the corrections would be insignificant at the level of precision of 
the data. 

The list of frequencies in Table 2 is pertinent to an interme- 
diate level of solar activity, which corresponds to a mean level 
of the 10.7-cm radio flux of 118 x 10~ 22 WnT 2 Hz -1 . The list in 
Table 3 is pertinent to the canonical "quiet-Sun" level of the ra- 
dio flux, which, from historical obse rvations of the index, is se t 
at 64 x 10~ 22 WnT 2 Hz ~ x (again, see lTapping & DeTracevlll99Cj) . 
To obtain low-degree frequencies pertinent to a level of activity 
higher than quiet-Sun, a linear interpolation can be performed be- 
tween the values in Table 2 and Table 3 (or a linear extrapola- 
tion may be carried out if the chosen activity level is higher than 
118 x 10~ 22 WirT 2 Hz 1 ). Both frequency tables have been pro- 
duced so that the frequencies correspond to the frequency centroids 
of the low-degree modes. The frequency centroids carry informa- 
tion on the spherically symmetric component of the internal struc- 
ture. 

The Sun-as-a-star data do not always provide estimates of the 
centroids by default. That is because only i + 1 of the 2^+1 
components in each non-radial mode are seen clearly in the ob- 
servations (those with even £ + m). At times of low solar activity, 
components in the non-radial modes are arranged symmetrically in 
frequency: then peak-bagging of the Sun-as-a-star data does give 
the required centroids. However, at times when there is significant 
surface activity - as in our 8640-d time series - the components 
are not arranged symmetrically in frequency, and so because some 
components are missing from the Sun-as-a-star observations we 
measure a slightly different frequency to the centroid. It is possi- 
ble to make a correction for th is effect, as discussed in detail by 
Appourchaux & Chaplinl Umk . We have applied such a correc- 
tion to give the 1 18 x 10~ 22 W irT 2 Hz 1 frequencies in Table 2. 

The frequencies in Table 3 are corrected to the quiet-Sun 
level, using a procedure that by definition corrects to the expected 
quiet-Sun centroid frequencies. The procedure assumes that the 
frequency correction can be parameterized as a linear function of 
the chosen activity proxy. The BiSON correction is self-consistent. 
First, we divide the long BiSON time series into shorter pieces in 
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Table 1. Raw, best-fitting frequencies (in /iHz) from our pseudo-global peak-bagging analysis. 



n 


£=0 


t=\ 


i= 


2 


t= 


3 




t= 


4 


6 


972.613 ± 0.002 


















7 
8 


1263.162 ±0.012 


1185.592 ±0.004 
1329.629 ± 0.004 


1394.680 


±0.011 












9 


1407.481 ±0.006 


1472.836 ± 0.005 


1535.861 


± 0.008 


1591.575 


± 


0.014 






10 


1548.336 ± 0.007 


1612.723 ± 0.007 


1674.540 


± 0.008 


1729.092 


± 


0.016 






11 


1686.601 ±0.011 


1749.285 ± 0.007 


1810.314 


± 0.009 


1865.307 


± 


0.019 






12 


1822.203 ±0.011 


1885.091 ±0.009 


1945.816 


±0.013 


2001.265 


± 


0.017 






13 


1957.431 ±0.012 


2020.810 ±0.010 


2082.131 


±0.015 


2137.821 


± 


0.019 






14 


2093.535 ± 0.013 


2156.812 ±0.014 


2217.698 


±0.018 


2273.563 


± 


0.026 


2324.163 


±0.051 


15 


2228.774 ± 0.014 


2292.032 ± 0.015 


2352.280 


± 0.017 


2407.707 


± 


0.025 


2458.716 


± 0.076 


16 


2362.823 ± 0.015 


2425.650 ± 0.014 


2485.948 


± 0.019 


2541.754 


± 


0.023 


2593.197 


± 0.053 


17 


2496.226 ± 0.016 


2559.235 ± 0.015 


2619.761 


±0.018 


2676.252 


± 


0.022 


2728.502 


± 0.038 


18 


2629.724 ± 0.014 


2693.437 ± 0.013 


2754.551 


± 0.016 


2811.440 


± 


0.020 


2864.349 


± 0.033 


19 


2764.211 ±0.014 


2828.258 ± 0.013 


2889.708 


± 0.016 


2947.075 


± 


0.018 


3000.227 


± 0.034 


20 


2899.101 ±0.012 


2963.447 ± 0.012 


3024.839 


± 0.016 


3082.439 


± 


0.022 


3136.105 


± 0.032 


21 


3033.845 ± 0.012 


3098.287 ± 0.013 


3159.997 


± 0.017 


3217.881 


± 


0.024 


3271.864 


± 0.048 


22 


3168.726 ± 0.014 


3233.321 ±0.016 


3295.275 


± 0.021 


3353.580 


± 


0.037 


3408.325 


± 0.069 


23 


3303.652 ± 0.019 


3368.689 ± 0.021 


3431.006 


± 0.030 


3489.699 


± 


0.049 






24 


3439.152 ±0.027 


3504.397 ± 0.028 


3567.248 


± 0.041 


3626.408 


± 


0.073 






25 


3575.083 ± 0.046 


3640.791 ± 0.036 


3703.388 


± 0.065 


3762.917 


± 


0.102 






26 


3710.938 ±0.086 


3777.436 ± 0.050 


3840.118 


±0.143 


3900.029 


± 


0.155 






27 


3847.250 ±0.177 


3913.995 ± 0.066 


3977.388 


± 0.297 












28 


3984.507 ± 0.323 



















Table 2. Corrected frequencies (in /iHz), pertinent to an intermediate level of solar activity (see text). 



n 




t= 


1 






t= 


2 






t= 


3 




6 


972.613 ±0.002 
























7 




1185.592 


± 


0.004 


















8 


1263.162 ±0.012 


1329.629 


± 


0.004 


1394. 


.680 


± 


0.011 










9 


1407.481 ±0.006 


1472.836 


± 


0.007 


1535 


.862 


+ 


0.007 


1591 


.575 


± 


0.014 


10 


1548.335 ± 0.007 


1612.723 


± 


0.007 


1674. 


.540 


± 


0.008 


1729 


.092 


± 


0.016 


11 


1686.601 ±0.011 


1749.285 


± 


0.007 


1810 


.314 


± 


0.009 


1865 


.306 


± 


0.019 


12 


1822.203 ±0.011 


1885.092 


± 


0.009 


1945 


.815 


± 


0.013 


2001 


.264 


± 


0.017 


13 


1957.433 ± 0.012 


2020.811 


± 


0.010 


2082. 


.130 


± 


0.015 


2137 


.820 


± 


0.019 


14 


2093.538 ± 0.013 


2156.813 


± 


0.014 


2217 


697 


± 


0.018 


2273 


.563 


± 


0.026 


15 


2228.779 ± 0.014 


2292.035 


± 


0.015 


2352. 


.280 


± 


0.017 


2407 


.709 


± 


0.025 


16 


2362.831 ±0.015 


2425.655 


± 


0.014 


2485. 


.950 


+ 


0.019 


2541 


.758 


± 


0.023 


17 


2496.237 ± 0.016 


2559.242 


± 


0.015 


2619. 


.766 


± 


0.018 


2676 


.259 


± 


0.022 


18 


2629.741 ± 0.014 


2693.450 


± 


0.013 


2754. 


.561 


± 


0.016 


2811 


.455 


± 


0.020 


19 


2764.233 ± 0.014 


2828.279 


± 


0.013 


2889. 


.726 


+ 


0.016 


2947 


.101 


± 


0.018 


20 


2899.133 ±0.012 


2963.478 


± 


0.012 


3024 


.867 


± 


0.016 


3082 


.471 


± 


0.022 


21 


3033.886 ± 0.012 


3098.327 


± 


0.013 


3160. 


.028 


± 


0.017 


3217 


.916 


± 


0.024 


22 


3168.773 ±0.014 


3233.358 


± 


0.016 


3295. 


.301 


± 


0.021 


3353 


.602 


± 


0.037 


23 


3303.698 ± 0.019 


3368.721 


± 


0.021 


3431. 


019 


± 


0.031 


3489 


.712 


± 


0.049 


24 


3439.195 ±0.027 


3504.417 


± 


0.028 


3567. 


.255 


± 


0.041 


3626 


.420 


± 


0.073 


25 


3575.121 ±0.046 


3640.799 


± 


0.036 


3703. 


.386 


± 


0.065 


3762 


.911 


± 


0.102 


26 


3710.972 ± 0.086 


3777.430 


± 


0.050 


3840 


.096 


+ 


0.143 


3900 


.013 


± 


0.155 


27 


3847.277 ±0.177 


3913.974 


± 


0.066 


3977 


.350 


± 


0.297 










28 


3984.527 ± 0.323 

























order to quantify the sensitivity of the frequencies to the changing SON observations was therefore also applied to the time series of 
10.7-cm flux (by linear regression). We then use this information the 10.7-cm radio flux, 
to calibrate our correction. Since we could not obtain robust es- 
timates of the t = 4 frequency shifts using our data, we did not 
attempt to correct the frequencies of these modes. There are gaps 
in the BiSON time series (e.g., due to inclement weather, and very 
occasionally instrument problems). The window function of the Bi- 

Full details of the correction procedures will be published 
elsewhere (Broomhall et al., in preparation). 
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Table 3. Corrected frequencies (in /iHz), pertinent to the canonical "quiet Sun" (see text). 



n £=0 t=\ t=2 i=3 



6 


972.613 ±0.002 








7 




1185.592 ±0.004 






8 


1263.162 ±0.012 


1329.629 ± 0.004 


1394.680 ±0.011 




9 


1407.481 ±0.012 


1472.835 ± 0.007 


1535.859 ±0.011 




10 


1548.333 ± 0.007 


1612.720 ± 0.007 


1674.534 ± 0.008 


1729.085 ± 0.016 


11 


1686.597 ±0.011 


1749.279 ± 0.007 


1810.304 ±0.009 


1865.293 ± 0.019 


12 


1822.196 ±0.011 


1885.081 ±0.009 


1945.797 ± 0.013 


2001.243 ±0.017 


13 


1957.421 ± 0.012 


2020.793 ± 0.010 


2082.102 ±0.015 


2137.788 ±0.019 


14 


2093.518 ±0.013 


2156.784 ±0.014 


2217.656 ±0.018 


2273.516 ±0.026 


15 


2228.749 ± 0.014 


2291.993 ±0.015 


2352.222 ± 0.017 


2407.643 ± 0.025 


16 


2362.788 ± 0.016 


2425.595 ± 0.015 


2485.873 ± 0.019 


2541.671 ±0.023 


17 


2496.180 ±0.017 


2559.162 ±0.015 


2619.668 ±0.018 


2676.149 ±0.023 


18 


2629.668 ± 0.015 


2693.347 ± 0.014 


2754.439 ± 0.017 


2811.318 ±0.021 


19 


2764.142 ±0.015 


2828.150 ±0.014 


2889.578 ±0.018 


2946.935 ± 0.021 


20 


2899.022 ± 0.013 


2963.322 ± 0.014 


3024.689 ±0.018 


3082.273 ± 0.025 


21 


3033.754 ± 0.014 


3098.140 ±0.015 


3159.821 ±0.020 


3217.683 ± 0.028 


22 


3168.618 ±0.017 


3233.139 ±0.018 


3295.063 ± 0.025 


3353.335 ± 0.041 


23 


3303.520 ± 0.021 


3368.469 ± 0.023 


3430.748 ± 0.034 


3489.409 ± 0.053 


24 


3438.992 ± 0.030 


3504.129 ±0.030 


3566.949 ± 0.044 


3626.078 ± 0.075 


25 


3574.893 ± 0.048 


3640.475 ± 0.038 


3703.044 ± 0.067 


3762.530 ±0.105 


26 


3710.717 ±0.088 


3777.067 ± 0.052 


3839.717 ±0.144 




27 


3846.993 ±0.177 


3913.570 ± 0.068 


3976.930 ± 0.298 




28 


3984.214 ±0.323 
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